/3-JACOBI PROCESSES 



N. DemnQ 

Abstract. We define and study a [0, l] m -valued process depending on three positive real 
parameters p,q,P that specializes for /3 = 1,2 to the eigenvalues process of the real and 
complex matrix Jacobi processes on the one hand and that has the distribution of the (3- 
Jacobi ensemble as stationary distribution on the other hand. We first prove that this process, 
called /3-Jacobi process, is the unique strong solution of the stochastic differential equation 
defining it provided that /3 > 0,p A q > m — 1 + 1//3. When specialized to /3 = 1,2, our 
results actually improve well known results on eigenvalues of matrix Jacobi processes. While 
proving the strong uniqueness, the generator of the /3-Jacobi process is mapped into the radial 
part of the Dunkl-Cherednik Laplacian associated with the non reduced root system of type 
BC. The transformed process is then valued in the principal Weyl alcove and this allows to 
define the Brownian motion in the Weyl alcove corresponding to all multiplicities equal one. 
Second, we determine, using stochastic calculus and a comparison theorem, the range of /3, p, q 
for which the m components of the /3-Jacobi process first collide, the smallest one reaches 
and the largest one reaches 1. This is equivalent to the first hitting time of the boundary of 
the principal Weyl alcove by the transformed process. Finally, we write down its semi group 
density. 

1. Introduction 

1.1. Motivation. The study of processes associated with root systems has known a considerable 
growth mainly after the birth of Dunkl operators which allowed to define Dunkl processes ([7]). 
The projection of the latter processes on Weyl chambers associated with reduced root systems, 
known as radial Dunkl processes, revealed close connections to eigenvalues of some matrix- valued 
processes ([!]) and to Brownian motions in Weyl chambers ([S],[I3], [IT], [IB])- Indeed, the 
eigenvalues of symmetric and Hermitian Brownian matrices are radial Dunkl processes associated 
with type A root systems while the square root of the eigenvalues of Wishart and Laguerre 
processes are radial Dunkl processes associated with type B root systems. Besides, in the complex 
Hermitian setting, the eigenvalues process is the l/-process in Doob's sense of a multidimensional 
Brownian motion killed when it first hits the boundary of the Weyl chamber (V stands for the 
Vandermonde function). Following the standard scheme in the theory of stochastic processes, 
one wonders how is the eigenvalues process of real and complex Jacobi matrix processes related 
to root systems and this was mainly behind our motivation. This connection is not new in its 
own and is implicitely written in [I] where authors relate the multivariate Gauss hypergeometric 
function defined via Jack polynomials to the hypergeometric function associated with the root 
system of type BC. Loosely speaking, the eigenoperator of the firstly-mentioned function is the 
generator of a diffusion that we shall define and call the /3-Jacobi process, where f3 > is the 
inverse of the Jack parameter (see [T] for more details and definitions). The name is justified 
by the fact that the 1-and 2 Jacobi processes fit the eigenvalues process of the real and complex 
Jacobi processes ([H]) and also from the fact that the stationary distribution of the /3-Jacobi 



Key words and phrases. Jacobi processes; strong uniqueness; principal Weyl alcove; root system of type BC; 
hittining time; multivariate Jacobi polynomials. 

^SFB 701, Fakultat fur Mathematik, universitat Bielefeld, Bielefeld, Germany, 
email: demni@math.uni-biclcfeld.de 



1 



process is given by the multivariate Beta distribution corresponding to the /?- Jacobi ensemble 
(0)- 



1.2. From matrix Jacobi processes to /3-Jacobi processes. The real (complex) matrix 
Jacobi process was introduced and studied in |llj . It was defined as the radial part of a left 
corner of a n x n orthogonal (unitary) Brownian motion. The latter, say O n , may be defined 
via the heat kernel in the orthogonal (unitary) group ([2]), and by a left corner of size m x p 
we mean the compression of 9„ by two projections P m ,Qp of ranks m < n,p < n, that is 
Prn®nQp- The real (complex) matrix Jacobi process of parameters (p,q = n — p) is then defined 
by J m '■= Pm®nQp®nPm where 0* stands for the adjoint of 0„. This definition is extended to 
real parameters p,q via stochastic differential equations. Let (3 = 1,2 respectively , then when 
P A q > (m — 1) + 2//3, the eigenvalues process is the unique strong solution of ([TT] p. 135): 
(1) 



d\i(t) = 2v/(Ai(£)(l - \i(t))dvi(t) + (3 



(p-(p + q)Mt)) + Y / 



A l (t)(l-A J (t)) + A J (i)(l-A l (t)) 



Ai(*)-M*) 



dt 



where (^i)™ i are independent Brownian motions and < A m < • • • < Ai < 1. In fact, the 
eigenvalues remain in ]0, 1[ when p A q > (m — 1) + 2/(3. When one considers a Haar orthogonal 
(unitary) matrix instead of an orthogonal (unitary) Brownian motion, the distribution of J no 
more depends on time (stationary) and its eigenvalues are distributed at any time t according to 
the multivariate Beta distribution: 



z^{\) : = c mf3 n \r 0(m - i)/2 (i ~ xtf-^-w n 

i—1 l<2<j<m 



for (3 = 1,2 respectively and for some normalizing constant C n , P ,q,p- In the same way the (3- 
Jacobi ensemble ([H]) generalizes the JOE, JUE, JSE ([BJ), we introduce the (3- Jacobi process 
as a solution, whenever it exists, of ([T]) for arbitrary (3 > 0. We will prove that (JlJ has a unique 
strong solution for p A g > (rri — 1) + 1/(3, for any time £ and for any starting point < A m (0) < 
• • ■ < Ai(0) < 1. This result actually improves the strong existence and uniqueness result derived 
in |llj . To proceed, we will instead focus on the process 4> defined by </>j := arcsin y/Xi, 1 < i < m 
so that < 4> m < ■ ■ ■ < 4>x < 7r/2. The process </> is shown to satisfy 

d(j)(t) = dv(t) +V V] fe(a) log(sin(o;, (f>(t)))dt = dv(t) + V] fc(a) cot ((a, (f>{t)))dt a, 



oefi 



+ 



where ^ = (Wj)^i, (■) is the usual inner product in ]R m , i? + is a positive system of the non 
reduced root system of type JSC and k : a h- > fc(a) is a positive multiplicity function, that is, 
k(a) = k(^) if and only if a, 7 are conjugated under the action of the reflections group (|15j). 
The strong uniqueness will follow from a scheme developed in [9] in order to prove a similar result 
for radial Dunkl processes. However, this setting needs more care since the state space of <j> (or 
A) is a bounded domain while Weyl chambers are not. In fact, 4>/tt is valued in the principal 
Weyl alcove A of type BC ([15J, p. 89) therefore one defines the Brownian motion in the Weyl 
alcove of type BC as the process ip/ir corresponding to k(a) = 1. Unlike Brownian motions in 
Weyl chambers and the Brownian motion in the Weyl alcove of type A (see below) , this process 
is not a /i-process of a m-dimensional Brownian motion killed when it first reaches the boundary 
dA of A , where 

h(4>) := ] J sin(< a, >). 

More precisely, h is shown to satisfy (A/2)h > ch, c < 0. The end of the paper is devoted to the 
semi group density of A expressed as a bilinear series of multivariate Jacobi polynomials ([20 ). 
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2. Strong uniqueness 



2.1. The transformed process. As mentioned before, the connection to root systems was 
implicitely written in [1] and was done via differential operators which are the generators of A 
and </>. For the reader's convenience, we display the details and derive the stochastic differential 
equation satisfied by <fi. Let < A m < • • • < Ai < 1 and consider fTJ) up to 



Let 



Using 



inf{ft, \ m (t) = 0} A inf{t, Ai(f) = 1} A {t, Xi(t) = A 3 ft)for some 
arcsin(VAi) := s(A,), 1 < i < m, then < m (O) < • ■ • < </>i(0) < ir/2 and: 

2 



s'(Xi 



1 



sin 2<fri 



2(2 sin 2 



1) 



2 cos 2cj)i 



sin J 2<fo 



sin 2 <f>i — sin 2 0j = sin(0i + 0j) sin(</>i — 4>j), 



sin 2 cos 2 d>j + cos 2 sin 2 6, — — fsin 2 ( 



fa) + sin 2 {fa- <£,•)]) 



then Ito's formula gives 
d</>ift) = <Mft) 

Using 



^ ( P -( P + g) S in 2 ^) _ cot2 ^ 



sin 20; 



, /3 eft 



E 



sin 2 (0, ft) + <^ft)) + sin 2 (0, ft) - <^ft)) 



1 — cos 20i 



sin 2(f>i 



2 sin2&(t) ^ sin(^ft) + ^(t)) sin(^ft) - fo(t)) 



cot ^ — cot 20j, 



then 



Using 



one gets 
1 



(p (p + g)sin 2 ^) _ = (pj cot(( ^. (t)) + ^ _ x] cot(2( ^. (t))> 

sin20i 2 



sin 20^ 
1 

sin 2 z 



[cot(0» + 0j) + cot(0j - sin(0i + </>j) sin(<^ - 0j), 

1 + COt 2 2, 



sin 2 (0 t ft) + 3 ft)) + sin 2 ft) - <j>j(t)) = l_[l/sin 2 (^ft) + &(<))] + [l/sin 2 (<^ft) - <^-ft))] 

2sin2<^(t) sin(^(t) + <£,•(*)) sin(<^ft) - <^-(t)) 2' 
_ 1 cot 2 (0, ft) + 4,j(t)) + cot 2 (0, ft) + + 2 



cot(^ft) + 0j-(t)) + cotOift) - ^(t)) 



2 cot ft) + 0,-(t)) + cot(^ft) - &(<)) 

1 - cot(^ft) + <j) 3 (t)) cot^ft) - <£,-(*)) COt(0ift) + <^ft)) + cot(^ft) - 3 -(t)) 



cot(^ft) + 4>j{t)) + cot(^ft) - ^(i)) 
Finally, using 

cot(w + u) 



cot (it) cot(u) — 1 
cot (it) + cot(t)) 



one gets 



where 
(3) 



l (t)=du l (t) 



fco cot (f>i + ki cot 2<f>i(t)dt + k,2 V^[cot(</>, + 4>j) + cot(^>j — </>j)] 



2k =P(p-q), ki = (3(q - (ro - 1)) - 1, 2fc 2 = /3. 
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Easy computations show that 7r/2 — <fi satisfies @ with (p, q) intertwined. 

2.2. Existence and uniqueness of a strong solution. From @, the generator of 4> ac ts 
on smooth functions as (it is the trigonometric version of the VF-invariant Dunkl-Cherednik 
Laplacian, [2"2"]): 

JSf :=iA-<V,V*(-)>, $(0):=- Y, ^(")logsin«a»), 

where 

R = {±e t , ±2et, 1 < i < to, ±(e* ± ej), 1 < i < j < to}, 
R+ = {e,, 2e l; 1 < « < m, (e» ± ej), 1 < i < j < to}, 
S = {ei-e i+ i,l<i<m-l,e m }, 

( e i)iLi is the canonical basis of K. m and S is the simple system corresponding to the positive 
system R + ([15]). R is known as the non reduced root system of type BC ([15] p. 41) and the 
action of the reflections group on R gives rise to three orbits so that the multiplicity function takes 
three values corresponding to {±ej, 1 < i < to}, {±2ei, 1 < i < to}, {±(ej±ej), 1 < i < j < to}. 
With regard to ([5]), fc(ej) = ho, 2k(2ei) = fci, k(ei±ej) = k,2- Setting fa := fa/ft, then is valued 
in the closure of the principal Weyl alcove ([15]. p. 90) defined by : 

A = {4> € R m , (a, fa > Ofor alia € S, (a, 0) < 1} = {4> E R m , < fa n < ■ ■ ■ < fa < 1/2} 

where a = 2ei is the unique highest root ([E] p. 40). When ko = p = q, (J2J) invovles the 
root system of type C given by the data ([H] p. 42) 

R = {±ei ± ej, 1 < i < j < m, ±2ej, 1 < i < to} 

= {ej ± ej, 1 < i < j < to, 2ei, 1 < i < to} 
S = {e, - e i+ i, 1 < i < to - 1, e TO } 

and we refer to the /3-Jacobi process as the /3-ultraspherical process. Now, we proceed to the 
proof of the strong uniqueness for ^ subject to k$,ki,k2 > 0. Once we did, we apply this 
result to 7r/2 — cf) and to the ultraspherical ,5-Jacobi process to deduce the strong uniqueness for 
p hq> (to — 1) + 1/(3. For, we need to briefly recall the scheme developed in ([3]) used to deal 
with radial Dunkl processes and inspired from [3J, [3]- First, a result from [3J, [1] states that 

d(j)(t) = dv{t) - V$(0(t)) + n{fat))dL t 

where n(fat)) is a unitary inward normal vector to A at fat) £ dA Q , that is ([3]) 

(4) {fat)~a,n{fat))) <0, VaeX, 

and L is the boundary process satisfying: 

dL t = l{0( t ) e9j4o }dL t , 

has a unique strong solution for all t > and 0(0) € A . Second, it remains to prove that L 
vanishes. For, we need to prove two Lemmas: 

Lemma 2.1. Let 

H a := {<P e R m , (a, fa — 0}, a G S, 
H & := {0eR m , (5,0) =^} 

6e i/ie walls of A . Then 

dA a = U aeSu{dl }H a n {fa (n(fa),a) ^ 0}, 
i/iat is if <p € dA , then there exists a £ S U {a} such that <f) S i? Q and (n((j)), a) ^ 0. 
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Once we prove this Lemma, the following trivially holds 

Lt = ^-{4>(t)adA }dL t < 2J l{<t>(t)&H a }1-{{n{<t>{t)),a)^o}dL t . 

a£Sll{a] 

As a matter of fact, the strong existence and uniqueness result will follow after proving that 
Lemma 2.2. 

1 {<P(t)eH a }(n{4>(t), dL t ) = 0, a£ Sl){a}. 

Proof of Lemma \2.1\ let <fr £ dA then either belongs to and one claims that n(<f>) = 
— ei = —a/2 is a unitary inward normal vector at <j>, therefore (n((f>),a) ^ 0. This claim follows 
immediately from ([U and from the fact that cpx = %/2. Or (<j) G 3j4 \ -ff<s) there only exists 
a £ 5 such that G and one proceeds as in Lemma 2 in [5] but with more care since the Weyl 
alcove is a bounded domain while the Weyl chamber is not. Hence, assume that for all a £ S 
such that (a, (j>) = 0, one has (n((f>), a) — 0. The idea is to find e > such that <t> — en(0) £ A 
and then conclude that n(<p) — by substituting a = <fi — en((f>) in 

If (n(<p), a) = for all a € S, then (n((/>), a) = and <f> — en(<fi) £ A for all e > and 4> £ dA . 
If (n(cp),a) ^ for some a £ S, then (a, 0) > 0. For those simple roots, either (a, n((j))) < 
thus (a,n(4>)) < (a £ R + ) and one seeks e such that 

{a,4>-en(<ff)) < 1. 

Such an e exists since we assumed that {&,</>) < 1. Finally, if for those simple roots, some are 
such that (a,n(4>)) > 0, then one seeks e such that 

(d>, a) 

< e < 



(n(<j)),a) 

for all a such that (0, a) > 0, {n(4>),a) > and such that —e(a,n((f>)) < 1 — (6t,4>). The first 
condition is satisfied provided that 

(th.rv) 

< e < 



(0,a)>O,(n(0),a}>O (n(0),a)' 

while for the second we need to write a = J2 a e 

> so that 

{a,n((j))) = ^2 a a (a,n(4>)) + a a (a,n(cp)) > ^ a a (a,n(4>)). 

<a,n(0))>O <Q,ri(0))<O (a,ra(0))<O 

Thus, 

-e(a,n(4>)) < -e 2J a a (a,n(4>)). 

(a,n(0)><O 

If there is no a such that (0, a) > 0, (n((j>), a) < 0, then (a, > since a £ i?+ and 

(a, — en(0)) < (a, </>) < 1. 

Otherwise 

0<e<[- ^ a a {a,n{cj)))]- l {l - (a,<f>)) 

(a,n(0))<O 

and one chooses the smallest e which satisfies both conditions. The Lemma is proved. 



Proof of Lemma \2.2\ when a £ S, the proof is exactly the same as in Lemma 1 in [9] when 
substituting x <— > — ln{a, x) by <p t— > ~ logsin(a, 0). Thus, one has to deal with the additional 

5 



term corresponding to Let 9((f>) := — logsin(</>), then the occupation density formula yields 



L a t (n ~ (a, <f>))\0 (a)\da = (a, a) / \6 (tt - (a, 0(a)))|da = 4 / |0 ((a, <£(s)))|ds 
o Jo Jo 

since cot(7r — z) = — cot(z), where L^(tt — (a, fa)) is the local time at time t and at level a of the 
real- valued semimartingale 7r — (a, fa). Then, following line by line the proof of Lemma 1 in [9], 
we get 

/>7T 

m-K- (a,fa))\6(a)\da < oo 

o 

which implies that L\ (tt — (a, fa)) = for all t since 0' blows up at 0. Then, we use Tanaka's 
formula to compute 

= d[ir - (a, <j>(t)) - (tt - (a, <l>(t))) + } = -l{n-a,<t>(t))=o}d(a, <f>(i)) 

and finally use similar arguments to prove the statement of the Lemma. H 
Thus we proved that 

Theorem 2.1. If k$ > 0, k\ > 0, &2 > 0, then (0) has a unique strong solution for allt > and 
0(0) G X- 

Since 7r/2 — cj> satisfies ([2]) with p, q intertwined and since the case p — q fco = involves 
the reduced root system of type C which share the same Weyl alcove and the same highest root 
with the root system of type BC, our proof applies to those cases. We finally deduce that 

Corollary 2.1. If p A q > (m — 1) + l/(3,/3 > 0, then (dp has a unique strong solution for all 
t > and (f)(0) G ~Al- 

Remarks. 1/Our result simplifies to p A q > m in the real case (3 = 1 and p A q > m — 1/2 in the 
complex one (3 = 2. 

2/The process <fi/ir may be interpreted as particles in the interval [0, 1/2]. For particles on the 
circle, see [5], [H]. In that case, <j) m < ■ ■■ < <p\ < 4> m + 2tt so that 4>/2ir is valued in the Weyl 
alcove of type A (the highest root is <\>\ — cp m , [15 ) p. 41). 

2.3. Brownian motion in Weyl alcoves. It is known that Brownian motions in Weyl cham- 
bers are radial Dunkl processes with a multiplicity function that equals one and that they are 
the unique strong solution of some stochastic differential equation with a rational singular drift 
for all t f[5].|13j). Similarly the Brownian motion in the principal Weyl alcove of type A or 
equivalently the eigenvalues process of a unitary Brownian motion satisfies (see the end of |14j ) 

d<f>i(t) = dvi(t) + cot(<fe(i) ~ <Pj{t))dt, 

which also has a unique strong solution for all time t (0). As a matter of fact, we define the 
Brownian motion in the principal Weyl alcove of type BC as the unique strong solution of (f2]) 
with k(a) = 1 for all a G R, that is ko = h-z = 1, k\ = 2 or (3 = 2, q = m+l/2,p = q+1 = m+3/2. 
Its generator reads 

^=\ A + E (V,Vlogsim>,0)) := + <V, V log /*(<£)>, 

a£R + 

where 

m 

h{4>) := Yl sm (( a ;0)) = J| sin(0 i ) sin(2^-) JJsin(^j - fa) sin(0 i + fa). 

aeB. + j=l j<l 
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Then, h is strictly positive on A and vanishes for (p € 0A o . Unfortunately, unlike the case 
of Brownian motions in Weyl chambers or the case of the Brownian motion in the principal 
Weyl alcove of type A, this process is not a /i-process of a mutidimensional Brownian motion. 
Indeed, we can prove (see Appendix) that h is not an eigenfunction of A/2 and rather satisfies 
Ah > ch for some strictly negative constant c. Nevertheless, when j3 — 2, A is known to be the 
V-transform (V stands for the Vandermonde function) in Doob's sense of m independent real 
Jacobi processes of special parameters conditioned never to collide p. 141). 

3. The first hitting time of dA 
We define the first hitting time of dA by 

T = mi{t > 0, {4>{t)/n) G dA a } = T a A inf T a 

where 

T a := inf{i > 0, < a, <j){t) > = 0}, 

T & := inf{* > 0, < a, <j>(t) >= 2<f> x = tt}, 

and (p is the unique strong solution for all £ > of with R — BC m and p /\q > (m — 1) + 1//3. 
We will prove via stochastic calculus that 

Proposition 3.1. 

• I/O < k 2 < 1/2 < j9 < 1, then T a < oo a.s. for a 6 {e, - e i+ i, 1 < i < m — 1}. 

• //0<fc + fci + l<2^0<p-(m-l)< 2//3, i/ien T em < oo a.s.. 

• < g — (to — 1) < 2/ ft, then = T 2ei < oo a.s.. 

Recall that for a radial Dunkl process X associated with a reduced root system R and a Weyl 
chamber C, we compared the one dimensional process (a, X), a G 5 to a Bessel process in order 
to prove the fmiteness of the first hitting time of dC ([H])- Here, we shall compare (a, (j>) , a G S 
to a real Jacobi process and use that ir/2 — <f> satisfies ([2]) to deal with T^. However, the case we 
have in hands is more delicate since the root system of type BC is non reduced. We start with 
recalling some needed facts about Jacobi processes from [TT], [21] . 

3.1. On real Jacobi processes. The real Jacobi process of parameters d,d! > is the unique 
strong solution of 

dJ t = 2y/,J t (l - J t )d lt + (d-(d + d')J t )dt, d, d' > 0, 

where 7 is a real Brownian motion. It is known from ( [11]) that J remains in ]0, 1[ when 
d A d' > 2, hits if < d < 2 and hits 1 if < 61 < 2. According to our result on strong 
uniqueness specialized to the rank one case BCi, "0 : — arcsin\/J is the unique strong solution of 

(5) dipt = d-yt + [(d - d') cot(^t) + {d! - 1) cat{2ip t )}dt. 

subject to d A d' > 1. 

3.2. Proof of the Proposition. Let ao G R+, then from (f2]), one has 

d < ao,(p(t) >— ||ao||c?7t + fc 2 1 |c^o 1 1 2 cot(ao, cf>(t))dt + \_. k(a)a(a) cot(a, (p(t))dt 

aER + \a 

where a(a) =< ao,a >. Now, let ao G {e^ — e<+i, 1 < i < m — 1} and let p > q. Denote o~q 
the reflection with respect to H ao . Then, for all a G R+ \ ao, easy computations show that 
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(cro(a), oq) = —a(a). Since <Jo(a) belongs to the same orbit containing a and is positive (see 
[TB] p. 10), one has 

F t :— k (a) a (a) cot (a, c/)(t)) = k (a) a (a) [cot (a, 4*{t)) — cot(<7o(a), <fi(t))]. 

a£R + \a a a£R + \ao 

a(a)>0 

As a result, one writes for any t > 0: 

d(<fc(t) - = v^dTi + 2* a cot(&(t) - cj) i+1 {t))dt + F t dt. 

This drift is strictly negative on {T ao =00} since <j) l— * cot</> is a decreasing function and since 
{a.Q,<f)(t)) > so that: 

(a - <ro(a), 4>{t)) = 2^(a , 0(i)> > 0. 
INI 

Using Proposition 2. 18. p. 293 and Exercice 2. 19. p. 294 in pj)] , one gets 

P m (Vt>0, (a ,<Kt)) <Z t ) = l 

where 

dZ t = V2d-y t + 2k 2 cot(Z t )dt, Z a = (a , 0(0)) 

is defined on the same probability space as <p. According to (0), (sin 2 Z t / 2 )t>o is a real Jacobi 
process of parameters d — 2k 2 + l,d' = 1 and the first statement of the Proposition is proved 
when p > q. Applying the same scheme to n/2 — <f> when q > p, we get the desired first statement. 
Now, let «o = e m . Compared to the previous case, one has to take {e m , 2e m } out of the drift 
term of {e m ,4>) to get that ao(a) £ R+ for a G R + \ {e m , 2e m }. The last fact is easily checked 
since for a = e, ± ej this amounts to consider the reduced root system of type B ([15] p. 10), 
otherwise e m is orthogonal to {ei,2ei, i 7^ to} so that a® fixes pointwise this set. Accordingly: 

d < ao, (f>it) >— d(f) m (t) — d-ft + k Q cot(<j) m {t))dt + k\ cot(2^ m (i)) + F t dt 

where 

F t = k(a)a(a)[cot((a,(f>(t)))-cot((a (a),<j>(t)))] 

a£R\\{e m ,2e m } 
a(a)>0 

where R\_ = {e^ — ej, 1 < i < j < to} since a(a) = for a € {ei,2ei, i 7^ m}, and fc(cc) = 0. 
Arguing as above, T Cm < 00 a.s. if0</c + &i + l<2<^0<p-(TO-l)< 2//3. 
Finally, recall that p := ir/2 — <j> satisfies © with (p, q) intertwined, that is 2fco = (3(q — p), fci = 
(3{p — (to — 1)) — 1,2^2 = P- Since /i m = 7r/2 — </>i and since, as shown above, the drift 
term Ft of d/j, m (t) — d(fj,(t),e m ) is independent from p, q, we conclude that p m hits a.s. if 
< j3(q - (m - 1)) < 2 thereby T 5 < 00 if < q - (to - 1) < 2/(3. U 

Remark. When translated from <fi to A, the previous statements specialize to results from |llj : 
in particular the eigenvalues process of both real and complex Jacobi processes never collide. 
However, it is known from [11] that the eigenvalues process never hit {0,1} when p A q > 
m — 1 + 2/(3 for (3 = 1, 2 respectively. This remains true for all (3 > by the same arguments 
used in [TT] , 
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4. Semi-group density 



We end this paper by writing down the semi group density of the /3- Jacobi process A. Before 
proceeding, we briefly consider two cases for which the semi group is known: the univariate case 
m = 1 and the multivariate complex Hermitian one corresponding to 2-Jacobi processes. In the 
former, the semi-group density reads (see [25] ) 



(6) 



Pt 



r(9,X) =^Te- 2r » t PZs(9)pr>°(\)W r > s (\), (M) € [0,1] 2 , r,s > -1, 



where p = 2(r+l),q = 2(a+l),r n = — n(n+r+s+l), (Pn' s ) n are orthonormalJacobi polynomials 
on [0,1] and W r ' s (X) := C r , s X r (l — X) s for some constant C r , s so that W defines a probability 
distribution. No closed form is known for this density, nonetheless an attempt to get a handier 
expression was tried in [TUj . In the latter, the 2- Jacobi process is the Vandermonde transform of 
m independent real Jacobi processes of parameters d — 2(p— m+1) := 2(r+l), d! = 2(q— m+1) := 
2(s + l) and conditioned never to collide p. 141). More precisely, the Vandermonde function 
V satisfies 



hV = — m(m— 1) 



2(m-2) d + d' 



2(m — 2) 

V = -m(m-l) ( - ! — ^ — '- +p + q- 2(m - 1) ] V := cV 



where L is the generator of m independent real Jacobi processes of parameters d, d! (see appendix 
in [TT] p. 219). It follows by Karlin-McGregor formula ([TH]) that the semi group density is given 

by 



p r t ' s '\e,X) :=e" 



V(9) 



dot [ ^z~ 2n{n+r+s+1)t K''\Q l ) p n s (^jW r ' s {x J ) 

W^{X) 



\n=0 



'det \Y,e~ 2n{n+r+s+1)t P^ s (9 l )P^ s (X j ) 



\n=0 



V{9)V{X) 



E E e 

&i £S m Til ,. . . ,n m >0 

E E 

cti ,(T2£5m n i > ■ ■ ■ > n m > 



2YhLi ni(m+r+s+l)t 



;.=1 



W^ S ' 2 (A) 



V(0)V(A) 

e -aEEa « oaW (». 3 co+H-.+i)t jj (ft)^ (A CTl(l) ) 



V(WA) 



where 



Tl" 



r,s,/3 



(A) := C,...,,,^ J J AI •: 1 - A,) s ]"[ |A< - A,f 



S m is the symmetric group, and {0 < A TO < • • • < Ai < 1}, {0 < 9 m < ■ ■ ■ < 9\ < 1}. Note that, 
for a given partition (n% > • • • > n m > 0) and f72 € S* m , one has 

m m 

E n g2(i)( n ff2(i) + r + s + 1) = E n *( n * + r + s+ 1). 

i=l i=l 

Thus summing first over a\ with the change of variables a = cri<72, one gets: 

2 ^,», f »,. +r+s+1H det[P^^-)]- =1 det[^(A 3 -)]^. =l T ^ a-2 

V{9) V(X) 



P r 2 (M) = e - c * 



ni>--->n m >0 



E 

— >n, 

E 



1)t det[^ s (^)]S=i det^^A,)] 



i .i= 1 urr,s,2 



ni>->n m >0 



7(9) 



y(A) 



Pt ,s ' 2 may be expressed as in © using the multivariate Jacobi polynomials P^ s ^ defined by 
Lassalle which take a determinantal form for (3 = 2 f |20p: 



p;< s < 2 (a) 



det[^_ i (A i )]« 



V(A) 

where r = (j\ > T2 > • • • r m ). Thus, performing the index change rii := + m — i, one finally 
gets 

(7) P? S ' 2 (M)= e-^P^(9)P^ 2 (X)W^ 2 (X), 



T 1 >...T m >0 



where we set 



rl 



n(n + r + s + 1 + 2(m - i)) 
»=i 

so that 

rn m 

Y ni(rii + r + s + 1) + c/2 = r 2 . + Y^(m - i)(m - i + r + s + 1) + c/2 
i=i i=i 

m— 1 m — 1 

= r 2 . + ^ i 2 + (p + g - 2m + 1) ^ i + c/2 

i=l i=l 

o (m — l)m(2m — 1) , m(m — 1) 

= r 2 T ~ K - >-+ > - + (p + q) 1 - ' +c/2 



With regard to © and |(ZJ), it is natural to claim that : 

Proposition 4.1. The semi group density of the (3- Jacobi process is given on {0 < X m < ■ ■ ■ < 
Ai < 1} by 

(8) p r t ' s ' P (0,\) ■= E e-^(J)^(W(A) 

ri>— >r m >0 

/or p A q > (to — 1) + l//3„ w/iere /3(p - (m - 1)) := 2(Y + 1), (3{q - (m - 1)) := 2(s + 1) and 

m 

:= E r^Tj + r + s + 1 + /3(m - i)). 

i=l 

Proof, given a bounded symmetric function / on [0, l] m , define for t > 

T t f(9):= f f(\)p r t ' s ' (9,\)d\ 

for e A' := {0 < #1 < ■ • • < m < 1} and T f = f. The above integral converges: this uses 
the boundness of /, the exponential term with strictly positive t and Fubini Theorem. Besides, 
Tfl = 1 and ||T t || is bounded for all t > 0. Indeed, for t > 0, the first claim follows from Fubini's 
Theorem, the orthogonality of (P^ S ^) T and Pq' s ' /3 = 1 so that the only non zero term is that 
corresponding to Tj = for all i. The proof of the second claim is almost similar. One also easily 
checks that T t T s = T t + S - Let Jzf# be the generator of the (3- Jacobi process 9. Then, P^ s '$ is an 
eigenfunction of Jzfg with eigenvalue 2r^ (|20J) so that J/fgT t f(9) = dtT t f{9) by the dominated 
convergence Theorem. Now, let us consider the Cauchy problem associated to ^£q: 

dtu f (t,9) =3?ou f (t,0) 

u f (0,-) = f, 
10 



where u f S C 1 ' 2 (R+ x K) n C b (JR + n X) with boundary condition: 

<Vu f (t,6),n(6) >=0 

where n(0) is a unitary inward normal vector at 8 <E <9A'. Define ut(f){0) :— Uf(t,9). It is 
known (|23p that the above Cauchy problem has a unique solution. Consequently u t — T t for 
alH > and (T t )t>o is the semi group of the /3-Jacobi process with density given by (JH]) (note 
that T t f — > / in L 2 -sense as i — * 0). ■ 



5. Appendix 

Below, we give a quick proof of the fact that Ah > ch for some strictly negative constant 
To proceed, recall that, for <p E A , 



H&) = JJ sin (0i) sin(2^i) JJsin(0i - <p 3 ) sin(<^ + <f>j) 

i—l i<j 

m 

= 2 m J]sin 2 (^)cos(^)n( sin2 (^) - sin2 (^)) = 2 m Y[X iy /l-X i V(X) 

i—l j<i i—l 

where A^ = sin for 1 < i < m, A = (A^)^ and V stands for the Vandermonde function. Set 

m 

«(A) = JTA. lV /l - X i: v(4>) = (sin 2 (<?!)i),- • • ,sin 2 (0 m )) = A 

i=l 

so that h{(j>) = 2 m {uV){X) = 2 m {uV)(v(cj))). It follows that 



and that 



Now, write 



dih(<t>) = 2 m [2sin(^)cos(^)]a i (uI/)(A) = 2 m [2^/X l (l - Xi)]di(uV)(X) 

d 2 M) = 2 m [2^X % {\ - Xi)]di[2y/Xi(l - Xijdi(uV)(X)}. 

= 2" l [4A l (l - Xi)df(uV)(X) + 2(1 - 2X i )d l (uV)(X)}. 



2 — 3A^ 2 — 3A^ 



2 V / T^A~ 2Aj(l-Ai) 



i(A) 



and ft 7 = 7ft log V so that 
2 - 3A, 



at(«v)(A) = 

and 

a? MO 



2A,-(1 - AO 



■ftlogF(A) 



MOW = 



A, 2(1 - Aj) 



ft log 7(A) 



(uV) (A) 



MO 



'(A) 



1 



1 

A7 " 2(1 -A,) 

1 1 



ft log 7(A) 



a 2 logI/(A)--2 



A 2 2(1 -A,) 1 



4(1 -A,) 2 A. t (l-A0 



+ [(ft log 7(A)) 2 + ft 2 log 7(A)] 



Thus 

-)2 



— W = { ~ 2 ^ i)(2 3Ai) + [2(1 - 2A0 + 4(2 - 3A0]ft logV(A) - 4 - + 4A,(1 - AO 



2-3A, 
Ai(l - AO 



A, 



ft log 7(A). 
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9 — 7\ 4- ^A 2 r) 2 T/ 
(10 - 16A0ft log 7(A) - 4 + At(1 _ A ^ + 4Ai(l - A^)^f(A). 



Noting that 2 - 1X { + 5X; = (Xi - l)(5Aj - 2) and using 



Ev"~* 1 \ -> A,- m(m— 1) 

di log V(X) = ^2 ~x~^x~l = ^o: = V 

4 m 

— K 4>) = -8m(m - 1) - 4m + £ A,(l - ,./,• > > - < - - ^ — ^ 



then 

= -8m(m - 1) - 4m + ^ £ Ml " h)dfV(X) - £ " 

Finally, from p. 219 in [TTj . 

2 £ A 4 (l - A i)9 ^(A) = !)(*»- 2) y(A) 

i=i d 

so that 

A/i.,, , , 4m(m— l)(m — 2) 2 9 4m(m— l)(m — 2) 

— (0) = -8m(m-l)-9m ^ -+£— >-8"~ J 

i=i 



m 2 — m : — — := c 



which finishes the proof. 

Remark. In [12] , the author comes to the 2- Jacobi process as a limit in distribution of a rescaled 
Markov chain on the so-called Gelfand-Tsetlin graph. 
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